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Abstract 



Constitutive equations are derived for the viscoelastoplastic response of amorphous glassy 
polymers at isothermal loading with small strains. A polymer is treated as an ensemble of 
cooperatively relaxing regions (CRR) which rearrange at random times as they are thermally 
agitated. Rearrangement of CRRs reflects the viscoelastic response of the bulk medium. At 
low stresses, CRRs are connected with each other, which implies that the macro-strain in 
a specimen coincides with micro-strains in individual relaxing regions. When the average 
stress exceeds some threshold level, links between CRRs break and relaxing domains begin 
to slide one with respect to another. Sliding of micro-domains is associated with the vis- 
coplastic behavior of polymers. Kinetic equations are proposed for viscoplastic strains and 
for the evolution of the threshold stress. These equations are validated by comparison with 
experimental data in tensile relaxation tests and in tests with constant strain rates. Fair 
agreement is demonstrated between results of numerical simulation and observations for a 
polyurethane resin and poly (methyl methacrylate). 

Key- words: Glassy polymers, Polyurethane, Poly(methyl methacrylate), Viscoelasticity, Vis- 
coplasticity 

1 Introduction 

A constitutive model is derived for the viscoelastoplastic response of amorphous glassy polymers 
at isothermal loading with small strains. The time-dependent behavior of glassy polymers has 
become the focus of attention in the past decade p|-pO|. This may be explained by the use 
of new experimental techniques which allow complicated time-dependent programs of loading 
to be carried out (loading and unloading with different strain rates interrupted by relaxation 
and superposition of small oscillatory strains on loading with constant strain rates) at various 
temperatures. 

The present study is concerned with the response of polymers in relaxation tests and in 
tests with constant strain rates. We confine ourselves to relatively slow programs of loading 
and to temperatures in the sub-T g region, where T g is the glass transition temperature. This 
implies that crazing and other phenomena associated with brittle fracture are excluded from 
our consideration, which focuses on the sub-yield and post-yield response of ductile polymers. 
To simplify the analysis, we deal with uniaxial loading at small strains. These limitations may 



be avoided by using conventional approaches, see, e.g., [21|, but we do not dwell on this issue. 
The purpose is to develop a constitutive model which describes the viscoelastic and viscoplastic 
response of amorphous polymers at the micro-level, on the one hand, and which is sufficiently 



simple to be used in the study of three-dimensional deformations of polymers with complicated 
geometry, on the other. We aim to design stress-strain relations which correctly predict two 
phenomena observed in experiments: (i) changes in elastic moduli and relaxation spectra induced 
by changes in the test's temperature and (ii) the influence of the strain rate and temperature 
on the apparent yield stress. 

To describe the viscoelastic response of glassy polymers, we employ the concept of cooperative 
relaxation [22]| . An amorphous polymer is treated as an ensemble of cooperatively rearranging 
regions (CRR). Any CRR is modeled as a globule consisting of scores of neighboring strands of 
long chains which change their position simultaneously because of large-angle reorientation of 
strands [23, 24 1. The characteristic length of a relaxing region amounts to several nanometers 
p5| . In the phase space, a CRR is treated as a point trapped in its potential well on the energy 
landscape. A CRR spends most time at the bottom level of its trap. At random times, it hops to 
higher energy levels being thermally agitated. We adopt the transition-state theory [26] which 
postulates that some liquid-like (reference) level exists on the energy landscape. When a CRR 
reaches this reference state in a hop, it totally relaxes. Viscoelasticity of a glassy polymer is 
modeled as sequential rearrangements of CRRs trapped in potential wells with various depths. 

CRRs are assumed to be connected with one another by some links (which model temporary 
crosslinks, entanglements and van der Waals forces). At relatively low stresses, the links prevent 
mutual displacements of CRRs, which implies that the micro-strains in relaxing regions coincide 
with the macro-strain in a specimen. When the average stress in an ensemble of CRRs exceeds 
some threshold value, links between CRRs break, and relaxing regions become (partially) free 
to change their positions with respect to each other. The mutual displacement of CRRs models 
the viscoplastic deformation at the micro-level. 

The exposition is organized as follows. Thermally induced rearrangement of CRRs is dis- 
cussed in Section 2. Section 3 deals with constitutive relations for the viscoelastic behavior of 
an amorphous polymer. Kinetic equations for the evolution of the threshold stress are proposed 
in Section 4. Section 5 is concerned with comparison between experimental data and results of 
numerical simulation. Some concluding remarks are formulated in Section 6. 



2 Kinetics of rearrangement 

The position of the liquid-like state at the glass transition temperature T g is accepted as the zero 
energy level on the energy landscape. The depth of a potential well where a CRR is trapped 
with respect to the zero energy level is denoted as w > 0. We suppose that the position of the 
reference state on the energy landscape may change with temperature and denote by fi(T) its 
ascent at some temperature T <T g with respect to the position of the reference level at T g (the 
function £l(T) satisfies the natural condition fi(T g ) = 0). 

Denote by Xq the concentration of traps per unit mass of the bulk medium. For simplicity, 
Xq is assumed to be a constant (independent of temperature and the loading history). Rear- 
rangement of relaxing regions is entirely characterized by the function X(t,r,w) which equals 
the concentration of CRRs (per unit mass at time t) trapped in cages with potential energy w 
which have last rearranged before instant r < t. The current distribution of traps is described 
by the probability density, p(t,w), of traps with energy w at time t. The functions p(t,w) and 
X(t, t, w) are connected by the formula 

X(t,t,w) 

P{t,w) = . (1) 

Xo 

In the sequel, we confine ourselves to equilibrated polymers with steady-state distribution of 
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traps. This means that the function p is independent of time. Bearing in mind that in the sub- 
T g region, the distribution of traps with various energies drastically depends on temperature T, 
we set p = p(T, w). 

Let q{uj)dw be the probability for a CRR to reach (in a hop) the energy level that exceeds 
the bottom level of its cage by some value u' located in the interval [u),u + duj\. Referring to 
1 27], we set 

q(u>) = aexp(— aw), 

where a is a material constant. The probability to reach the reference state in an arbitrary hop 
for a polymer equilibrated at a temperature T reads 



Q(7» 



w+n(T) 



q{bj)doj = exp 



-a(w + Q(T) 



Denote by Tq = Tq(T) the attempt rate (the average number of hops in a trap per unit time). 
The rate of rearrangement R equals the product of the attempt rate Tq by the probability of 
reaching the reference state in a hop Q, 



R{T,w) = r(T) exp(— aw) 
where the relaxation rate T(T) is given by 

r(T) =r (T)exp[-aft(T) 



(2) 



(3) 



Equating the relative rates of rearrangement to the function R, we arrive at the differential 
equations 



1 



dX 



X(t,0,w) dt 



(t,0,w) = -R(T,w), 



dX 
~d7 



1 -1 ^2 



(t,T,w) 



& 2 x 

dtdr 



(t,T,Ul) 



-R(T,w). (4) 



X(t,0,w) 

dX 
dr 



The solutions of Eq. (4) with the initial condition X(0,Q,w) = Xqp(T,w) [see Eq. (1)] read 

Xqp(T, w) exp —R(T,w)t 
-(t,T,w) = XqF(t, w) exp —R(T,w)(t — T) 

where 



(5) 



. . 1 dX . 

F ^ w) = x- ^ T ^ 



t=T 



The number N(t, w) of relaxing regions (per unit mass) in traps with potential energy w rear- 
ranged within the interval [t, t + dt] equals the sum of the numbers of CRRs rearranged for the 
first time at instant t and those that have last rearranged at some instant r G [0, t) and reach 
the reference state at time t, 



N(t , w) = _» (t , ,„) -J 



t d 2 x 

dtdr 



(t, r, w)dr. 



(6) 



Neglecting the duration of a hop (a few picoseconds [23]) compared to the characteristic time 
of relaxation in the sub-T g region, we find that N(t, w) coincides with the number of CRRs 
landing within the interval [t, t + dt] in traps with potential energy w, 



N(t,w) = X F(t,w) 



(7) 
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Substitution of Eqs. (4) to (6) into Eq. (7) results in the Volterra equation 



F(t,w) = R(T,w)^p(T,w)exp -R(T,w)t + J F(r,w)exp -R(T,w)(t- 
The solution of Eq. (8) is given by 

F(t,w) = R(T,w)p(T,w), 
which, together with Eq. (5) implies that 
dX 

— (t, t, w) = X R(T, w)p(T, w) exp [-R(T, w)(t - r) 

3 Constitutive equations 



dr 



(8) 



(9) 



We assume that the macro-strain e equals the sum of the viscoelastic strain e c and the viscoplastic 
strain e p , 

e(t)=ee(t)+e p (t). (10) 

The viscoelastic strain, e e , is the average strain in relaxing regions, and the viscoplastic strain, e p , 
is the average strain induced by displacement of CRRs one with respect to another. To simplify 
our analysis, we neglect the distribution of micro-strains and suppose that the micro-strain in 
any CRR coincides with e e . 

A CRR is assumed to totally relax when it reaches the liquid-like state, which means that 
the natural (stress-free) configuration of a CRR after rearrangement coincides with the actual 
configuration of the bulk medium at the instant of rearrangement. For uniaxial deformation, 
the strain from the natural configuration of a relaxing region to its actual configuration at time 
t is given by 

e(t,r) = ee(t)-€ e (r), (11) 

where r < t is the last instant when the region was rearranged. A CRR is modeled as a linear 
elastic medium with the mechanical energy 

<K*> r ) = 2 c£2 (*' r )' 

where c = c(T) is the rigidity of a relaxing region. To minimize the number of material param- 
eters in the constitutive equations, we assume the rigidity c to be constant and independent of 
the energy of trap w. Summing the mechanical energies of CRRs and neglecting the energy of 
their interaction, we find the strain energy density of a polymer (per unit mass) 



*(t) = \c 



! (t,0) / X(t,0,w)dw+ e 2 (t,r)dT —{t,T,w)dw 

JO JO Jo OT 



Substitution of expressions (10) and (11) into this equality results in the formula 
1 



*(t) 



(e(t) - e v^)f J Q X(t,0,w)dw 



dX 



(t, r, w)dw 



(12) 
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At small strains, the average stress a(t) is expressed in terms of the macro-strain e(t) by the 
conventional formula 

- »W- (13> 

where p is mass density in the stress-free state. Combining Eqs. (12) and (13), we arrive at the 
stress-strain relation 



a(t) = Eo(T) 
+ 



(e(t) - e p (t)) y X(t,Q,w)dw 
o ((e(*)-e(T))-(e p (t)-e p (r)))dT jf ■^(i.T, 



(14) 



where Eq{T) = pc(T)Xo is the initial Young modulus. 

To proceed with the analysis of Eq. (14), a concrete form of the distribution function 
p(T,w) should be established. Referring to the random energy model [pq| , we assume that the 
probability density p is quasi-Gaussian, 



p(T, w) = C (T) exp 



(w - W(T)) 



2l 



2£ 2 (T) 



(to>0), P(T,w)=0 (w<0), (15) 



where W{T) and S(T) are parameters of the Gaussian distribution and the constant C is found 
from the condition 

poo 

I p(T, w)dw = 1. 
Jo 

4 Kinetics of plastic flow 

To describe the evolution of the plastic strain e p (t), we propose constitutive equations similar 
(to some extent) to the viscoplasticity theory based on overstress (back stress pQfl ), see 
also recent publications []17], [OJ. For simplicity, we confine ourselves to active loading with 
e(t) = de{t)/dt > 0, which means that the internal parameter A (the material time associated 
with loading) may be set equal to the macro-strain e. Our approach is based on the following 
hypotheses: 

1. The level of plastic strain is determined by the threshold stress g which depends on tem- 
perature and the loading history. 

2. Plastic deformation occurs only when the current average stress a(t) exceeds the threshold 
stress g(t). The rate of change in the viscoplastic strain e p (with respect to the internal 
time A) is proportional to the difference a — g, 

^L = Kl H(a-g), (16) 



where K\ is a material parameter and TL{x) is the Heaviside function, 

H(x) = 



1, x>0, 
0, x < 0. 



3. For a stress-free material, the threshold stress has a fixed value go. At the micro-level, go 
determines the strength of links between CRRs which prevent their mutual displacements. 
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In the general case, the quantity go depends on temperature T and on the strain rate e 
with which a specimen is loaded. A decrease in go with temperature is explained within 
the framework of the theory of thermally activated processes: the growth of temperature 
implies an increase in the amplitude of thermal fluctuations, which reduces the strength 
of links between CRRs. The effect of strain rate is also quite natural provided that the 
links break according to the brittle pattern. 

4. When the current stress a exceeds the threshold stress, the quantity g changes. The rate 
of change in the threshold stress is determined by two mechanisms. The first is associated 
with a decrease in the current value of g caused by partial destruction of links between 
CRRs (we suppose that some distribution of the links' strength exists, and the links with 
low strengths break first). The rate of decrease in the threshold stress (with respect to the 
internal time A) is assumed to be proportional to the difference a — g. Bearing in mind 
that g is a nonnegative function, we set 

^- = -K 2 gH(a-g), (17) 
where K 2 is a constant. Equation (17) ensures that g(t) = for any t > to, provided that 

g(t ) = o. 

5. The other mechanism for changes in the threshold stress g is associated with thermally 
induced healing of links between CRRs. We suppose that this process is similar to craze 
healing in glassy polymers, see |yj and the references therein. Because the healing process 
is driven by thermal fluctuations, its kinetics is analogous to that for the reformation 



process in a temporary network [32]. This means that the rate of reformation (per unit 

tional to the number of broken links, 

K 3 H(g -g), (18) 



time t) for links between CRRs is proportional to the number of broken links, 

dg 



dt 

where K 3 is a constant. 

6. In general, the parameters (i = 1,2,3) depend on temperature T and the loading 
history. To reduce the number of adjustable parameters in the model, we assume that the 
effect of these factors is rather weak and treat JQ as material constants. 

When A coincides with the macro-strain e, Eqs. (16) to (18) are transformed into the nonlinear 
differential equations 

^ = K x H{a - g)±, | = K 3 H(g - g) - K 2 gH(a - *)| (19) 
with the initial conditions 

e p (0) = 0, g(0) = g . (20) 

Given a loading program, e(t), Eqs. (2), (14), (15), (19) and (20) uniquely determine the average 
stress a(t). The constitutive equations are determined by four functions of temperature, Eq(T), 
W(T), S(T) and T(T), one function of temperature and strain rate, go(T, e), and three constants 
K\, K 2 and K% (without loss of generality, the parameter a in Eq. (2) may be set to unity). The 
total number of adjustable parameters in the model is essentially less than that in conventional 
stress-strain relations in viscoelastoplasticity based on the overstress concept Q 17, 18 1. 
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5 Comparison with experimental data 



To verify the constitutive equations, we determine adjustable parameters by fitting observations 
in tensile relaxation tests and in tensile tests with constant strain rates. We begin with the 
standard isothermal relaxation test with the loading program 




where eo is a given strain (which is assumed not to exceed the yield strain). It follows from Eqs. 
(19) and (20) that 

e p (t) = 0, g{t) = g {T,0). 
Substituting these expressions and Eqs. (2), (15) and (21) into Eq. (14), we find that 

■(w- W(T)) 2 



POO 

E(T,t) = C(T)E (T) exp 
J o 



2£ 2 (T) 



+ T(T)exp(-w)t 



dw, 



where 

E(T,t) = ^ 
eo 

is the current Young modulus. 

We begin with matching experimental data for a polyurethane resin (T g = 80 °C). A descrip- 
tion of specimens and the experimental procedure can be found in [p^ j. First, we fit observations 
in relaxation tests at various temperatures T. Given a parameter Eq, the quantities W, £ and 
r are found by using the steepest-descent algorithm. Afterwards, the initial Young modulus Eq 
is determined by the least-squares technique. 

Figure 1 demonstrates excellent agreement between experimental data and results of numer- 
ical simulation. The parameters W, S, Eq and T are plotted versus the degree of supercooling 
AT = T g — T is Figures 2 and 3. These figures show that far below the glass transition temper- 
ature, the quantities W, X, Eq and T are fairly well approximated by the "linear" functions 

W = a + a x AT, E = b + &iAT, E = c + a AT, log T = d + d\AT (22) 

with adjustable parameters aj, hi, Ci and d{. The values of these quantities at T = 75 °C, i.e. 
in the close vicinity of the glass transition temperature, deviate from dependences (22). This 
may be explained by the fact that at a very small degree of supercooling (AT = 5 K), the 
relaxation test was carried out at the strain (e = 0.01) that exceeded the yield strain. As a 
result, experimental data (curve 4 in Figure 1) demonstrate a combined effect of viscoelastic 
and viscoplastic deformations. This issue has been previously discussed in @. 

The parameters W and £ increase with a decrease in temperature, which is in good agree- 
ment with an assumption about an increase in the ruggedness of the energy landscape with 



the growth of the degree of supercooling [34]. The initial Young modulus Eq decreases with 
temperature in accord with conventional observations for the effect of temperature on elastic 
moduli. Surprisingly, the relaxation rate T decreases with temperature. This conclusion does 
not, however, contradict the theory of thermally activated processes, because the relaxation rate 



T is defined [see Eq. (3) 
ature) by the factor exp 



as the product of the attempt rate Tq (which increases with temper- 



-O(T) 



that characterizes the ascent of the reference energy level at 



temperature T with respect to that at the glass transition temperature. 
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We now approximate experimental data in a tensile test with a constant strain rate eo, 




<*) = ^ n (23) 



For this purpose, we integrate Eqs. (14), (15), (19) and (20) with the material parameters 
W, £ and T determined in relaxation tests and find the parameters go and Ki which minimize 
deviations between results of numerical simulation and experimental data. Because the initial 
Young modulus -Eo found in relaxation tests underestimates stresses at the initial region of 
the stress-strain curve, the value of Eq was found by fitting experimental data in the interval 
e G [0,0.02]. 

First, we approximate the stress-strain curve at room temperature (T = 25 °C) obtained at 
the strain rate eo = 0.01 min -1 . The quantities Ki and go are found by applying the steepest- 
descent procedure. Afterwards, we fix the parameters Ki and match experimental data at other 
temperatures and other strain rates with the only adjustable parameter go. Figures 4 and 5 
demonstrate good agreement between experimental data and results of numerical analysis. The 
initial threshold stress go is plotted versus the strain rate eo in Figure 6. This figure reveals 
that the dependence go(e) is fairly well approximated by the function (conventionally used to 
describe the effect of strain rates on the yield stress) 

go = /3o + /?iloge (24) 

with adjustable parameters 

We now repeat fitting experimental data for poly(methyl methacrylate) (T g = 105 °C) at 
room temperature. Because of the lack of observations in relaxation tests and in tests with con- 
stant strain rates for the same specimens, we use experimental data in relaxation tests exposed 
in |35| (the relaxation master curve corresponding to various aging times after quench from 
165 °C) and experimental data in tests with constant strain rates measured in We begin 
with matching observations in tensile relaxation test (21) with the strain e = 0.005. Figure 7 
demonstrates fair agreement between experimental data and results of numerical simulation. 
Using material parameters found in the relaxation test, we approximate observations in tests 
with constant strain rates. By analogy with the analysis of experimental data for polyurethane, 
we determine the coefficients Ki by fitting data in a test with the strain rate eo = 0.01 s _1 , fix 
their values, and match observations for other strain rates by using the only adjustable parame- 
ter go- Figure 8 shows excellent agreement between experimental data and results of numerical 
analysis. The dependence go(i) is depicted in Figure 6, which reveals the high level of accuracy 
of Eq. (24). 

To assess the effect of temperature on the viscoelastoplastic response of poly(methyl metha- 
crylate), we fit observations in relaxation tests (data are adopted from |37f ) and in tests with a 
constant strain rate (measurements are taken from [^]) at various temperatures in the sub-T g 
region. Figure 9 reveals good agreement between experimental data and results of numerical 
simulation. Because experimental data are available in a relatively small interval of time (3 
decades), we reduce the number of adjustable parameters by setting W = (which implies 
that only three constants are employed to match observations). The functions E(T) and T(T) 
are plotted in Figure 10, which demonstrates that the parameter £ increases with the degree 
of supercooling (in agreement with results presented in Figure 2). The rate of relaxation T 
decreases with AT, in agreement with the theory of thermally activated processes. It is worth 
noting that the parameter £ vanishes (the energy landscape becomes homogeneous) at the 
critical temperature T cr = T„ + 4.7 K, which is quite close to the critical temperatures for other 
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amorphous polymers, see 38]. The initial Young modulus Eq is plotted versus temperature T 
in Figure 11. The dependence Eq(AT) is correctly described by Eq. (22), where the coefficient 
Cj take different values in the close vicinity of T g , AT < 20 K, and far below the glass transition 
point. Observations in tests with constant strain rates are depicted in Figure 12 together with 
results of numerical simulation (the initial Young modulus Eq(T) was taken directly from the 
data in relaxation tests). This figure also reveals fair agreement between experimental data 
and the model's predictions. Material parameters go and K{ for polyurethane and poly(methyl 
methacrylate) are listed in Table 1 which shows that these quantities have the same order of 
magnitude for the two polymers. 



6 Concluding remarks 

Constitutive equations have been derived for uniaxial isothermal response of glassy polymers 
at small strains. An amorphous polymer is treated as an ensemble of CRRs trapped in their 
cages. At low stress levels, relaxing regions are connected one with another, which ensures that 
macro-strain in a specimen coincides with the average micro-strain in a CRR. When the average 
stress exceeds some threshold level, links between CRRs break, and relaxing regions can slide 
one with respect to another (which induced viscoplastic deformation at the micro- level) . Simple 
kinetic equations are proposed for changes in the threshold stress and plastic strain. The stress- 
strain relations are verified using experimental data for a polyurethane resin and poly(methyl 
methacrylate) at various temperatures and strain rates. Fair agreement is established between 
experimental data and results of numerical simulation. 
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Figure legends 



The Young modulus E GPa versus time t s for polyurethane at temperature T °C. Circles: 



experimental data [33]. Solid lines: results of numerical simulation. Curve 1: T = 25.0; 



curve 2: T = 35.0; curve 3: T = 50.0; curve 4: T = 75.0 

2. The dimensionless parameters W (unfilled circles) and £ (filled circles) versus the degree 
of supercooling AT °C for polyurethane. Circles: treatment of observations [03]] . Solid 
lines: approximation of the experimental data by Eq. (22). Curve 1: ao = —12.8947, 
d = 0.7053; curve 2: b = -9.4737, b x = 0.4263 

3. The initial Young modulus Eq GPa (unfilled circles) and the relaxation rate T s _1 (filled 
circles) versus the degree of supercooling AT °C for polyurethane. Circles: treatment 
of observations [33]. Solid lines: approximation of the experimental data by Eq. (22). 



Curve 1: c = 1.8221, c x = 0.0100; curve 2: d = -1.9310, d x = 0.0571 
4. Stress a MPa versus strain e for polyurethane loaded at temperature T °C with the strain 



rate eo = 0.01 min 1 . Circles: experimental data [33]. Solid lines: results of numerical 
simulation. Curve 1: T = 25.0; curve 2: T = 50.0 

5. Stress a MPa versus strain e for polyurethane loaded at T = 25 °C with the strain rate 
eo min -1 . Circles: experimental data |]33] 1. Solid lines: results of numerical simulation. 
Curve 1: eo = 0.1; curve 2: eo = 0.01; curve 3: eo = 0.005 

6. The initial threshold stress go versus the rate of strain eo min" 1 at T = 25 °C. Circles: 
treatment of observations for a polyurethane resin |33| and poly(methyl methacrylate) [|3G|]. 
Solid line: approximation of the experimental data by Eq. (24). Curve 1: PU, ao = 88.096 
and a% = 13.448; curve 2: PMMA, a = 97.922 and a x = 20.176 

7. The Young modulus E GPa versus time t s for poly(methyl methacrylate) at T = 25 °C. 
Circles: experimental data |35|]. Solid line: results of numerical simulation with Wo = 12 A, 
S = 14.8, E = 2.0442 GPa and T = 0.05 s -1 

8. Stress a MPa versus strain e for poly(methyl methacrylate) loaded at T = 25 °C with 
the strain rate eo s _1 . Circles: experimental data |36|]. Solid lines: results of numerical 
simulation. Curve 1: eo = 10 -2 ; curve 2: eo = 5 x 10~ 3 ; curve 3: eo = 10~ 3 ; curve 4: 
eo = 5 x 10~ 4 ; curve 5: eo = 10 -4 

9. The Young modulus E GPa versus time t s for poly(methyl methacrylate) at temperature 



T °C. Circles: experimental data [37]. Solid lines: predictions of the model. Curve 1: 



T = 50; curve 2: T = 75; curve 3: T = 85; curve 4: T = 90; curve 5: T = 100 

10. The dimensionless parameter £ (unfilled circles) and the rate of relaxation T s" 1 (filled 
circles) versus the degree of supercooling T °C for poly(methyl methacrylate). Circles: 
treatment of observations p7[ . Solid lines: approximation of the experimental data by Eq. 
(22). Curve 1: b = 3.4931, b x = 0.7483; curve 2: d = 1.0853, d x = -0.0564 

11. The initial Young modulus Eq GPa versus the degree of supercooling AT °C for poly(meth- 



yl methacrylate) . Circles: treatment of observations [37]. Solid lines: approximation of the 
experimental data by Eq. (22). Curve 1: Co = 0.3217, c\ = 0.0504; curve 2: Co = 1.1763, 
d = 0.0069 
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Stress a MPa versus strain e for poly (methyl methacrylate) loaded at temperature T °C 

-l 



with the strain rate eo = 0.1 min 
of numerical simulation. Curve 1: 



Circles: experimental data [16]. Solid lines: results 
T = 50; curve 2: T = 75 
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Table 1: Adjustable parameters go and Ki for polyurethane and poly(methyl methacrylate) 



Material 


T °C 


90 MPa 


Ki MPa" 1 


K 2 MPa" 1 


K% min 1 


PU 


50.0 


64.0 


0.02 


0.8 


35.0 


PU 


75.0 


38.0 


0.02 


0.8 


35.0 


PMMA 


50.0 


36.0 


0.07 


1.0 


20.0 


PMMA 


75.0 


32.5 


0.07 


1.0 


20.0 
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